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摘 要 利用 穷 举 搜索 法 结合 参数 估计 方法 确定 地 下 岩 土 热 物性 参数 ， 通 过 模拟 验证 该 方法 的 可 靠 性 ,并 用 于 实际 工程 
测量 ， 且 与 优化 方法 计算 结果 进行 对 比分 析 . 结果 表明 穷 举 搜索 法 可 以 在 搜索 区 域 得 到 目标 函数 方差 和 最 小 值 , 基于 此 得 
到 热 物 性 参数 值 可 靠 性 较 高 ， 而 优化 方法 得 到 的 目标 函数 方差 和 为 局 部 区 域 的 极 小 值 ， 其 确定 的 土壤 导热 系数 可 靠 性 较 


高 , 但 钻 孔 内 热 阻 和 容积 比 热 计算 结 果 可 信 度 相对 较 低 . 
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埋 管 中 循 环流 动 ， 把 地 下 浅 层 土壤 作为 热源 ,并用 
其 温度 随 季节 变化 浮动 较 小 的 特点 ， 将 循环 水 作为 
载 热 流体 可 在 冬季 从 其 提取 热量 、 夏 季 从 其 提取 准 
量 供给 用 户 四 ,地 源 热 泵 地 埋 管 的 换 热能 力 与 土壤 
的 热 物性 参数 有 着 直接 关系 ，Kavanaugh 的 研究 显 

中 ， 当 地 下 土壤 的 导热 系数 或 热 扩 散 率 发 生 10% 
的 偏差 时 ， 地 下 埋 管 设计 长 度 偏差 为 45%~5.8%， 
从 而 影响 初 投资 。 因 此 准确 确定 土壤 热 物 性 参数 是 
合理 设计 经管 孔 深 和 个 数 的 前 提 1 

目前 便 热 流 测试 法 四 是 国际 上 普遍 采用 的 获 
取 地 下 土壤 热 物性 参数 值 的 方法 ， 利 用 该 方法 结合 
建立 的 地 埋 管 与 周围 土壤 传 热 模型 ， 确 定 地 下 土壤 
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的 导热 系数 、 容 积 比 热 容 等 热 物 性 参数 中 .常用 的 
传 热 模型 有 线 热源 模型 四、 空心 柱 热源 模型 四 
心 柱 热源 模型 00 等 ， 计 算 土 壤 热 物性 参数 值 通常 
采用 斜率 法 UU 和 参数 估计 法 M, 斜率 法 可 以 方便 
地 得 到 土壤 导热 系数 ,但 不 能 确定 容积 比热容 1, 
参数 估计 法 可 以 同时 确定 导热 系数 、 容 积 比 热 以 及 
ATL PI TRECE, 在 工程 计算 中 被 广泛 应 用 04719), 
参数 估计 法 通常 是 利用 模型 计算 得 到 循环 水 温度 并 
将 其 与 实验 测量 结果 对 比 ， 寻 找 不 同 热 物性 参数 组 
合 值 条 件 下 两 者 之 间 的 方差 和 最 小 值 ， 方差 和 取得 
最 小 值 时 所 对 应 的 土壤 热 物性 数值 即 认为 是 测试 得 
到 的 准确 热 物性 参数 值 . 

利用 传 热 模型 公式 结合 参数 估计 法 求解 热 物 性 
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参数 的 过 程 实际 上 是 求解 多 元 非 线性 方程 的 过 程 ， 
许多 数学 家 致力 于 利用 优化 方法 更 迅速 地 找到 符合 
要 求 的 解 ,， 芝 用 的 传统 方法 主要 有 和 牛顿 法 、 弦 割 法 、 
抛物 线 法 等 ， 但 这 类 方法 需要 假定 初始 参数 值 ， 且 
最 终结 果 受 初始 参数 值 选 定 的 影响 很 大 . 遗传 算法 、 
蚁 群 算 法 等 避免 了 结果 对 初始 值 的 强烈 依赖 性 ， 虽 
然 寻 优 速度 较 快 , 但 并 未 实现 真正 意义 上 的 全 局 范 
围 内 的 搜索 ， 所 以 很 可 能 将 某 局 部 区 域 的 极 小 值 作 
为 全 局 范围 内 的 最 小 值 输出 91. 本 文 尝试 采用 穷 举 
搜索 法 U, 在 土壤 热 物性 参数 值 可 能 的 范围 内 寻找 
目标 函数 最 小 值 , 实现 全 局 范围 内 的 搜索 , 进而 确定 
热 物 性 参数 ,并 与 优化 方法 的 结果 进行 对 比分 析 . 


1 地 埋 管 与 周围 土壤 传 热 模型 
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出 口 位 置 循环 水 温度 随时 间 的 变化 ， 得 到 从 测试 开 
始 到 热平衡 过 程 中 温度 随时 间 的 响应 曲线 .根据 测 
试 得 到 的 数据 结合 地 埋 管 与 周围 土壤 传 热 模 型 进行 
计算 ， 并 利用 参数 估计 法 最 终 得 到 热 物性 参数 值 是 
目前 普遍 的 热 物 性 参数 计算 方法 ， 

由 于 工程 中 普遍 采用 的 线 热源 模型 计算 结果 满 
足 准 确 性 要 求 ， 且 较 其 他 模型 计算 便捷 、 适 用 人 性 
强 gj， 本 文采 用 线 热源 模型 作为 地 埋 管 与 地 下 周 
围 土壤 换 热 的 计算 模型 。 将 地 埋 管 看 做 线 热 源 与 外 
界 换 热 需 做 如 下 假设 P: 1) 地 下 土壤 视 为 无 限 大 介 
质 ; 2) 地 埋 管 位 于 钻 孔 中 心 位 置 ， 并 视 为 一 根 无 限 
长 的 线 热源 ; 3) 钻 孔 内 回填 材料 及 地 下 周围 土壤 各 
向 同性 ， 且 物性 参数 维持 恒定 ; 4) 锁 孔 内 不 存在 热 
基 累 积 现象 ; 5) 忽略 土壤 内 的 水 分 迁移 、 不 同 材质 
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恒 热 流 测试 法 是 在 现场 进行 热 响应 实验 ,采用 R s 
I 恒定 功率 加 热 地 埋 管 循环 水 , 循环 水 在 埋 管 中 循环 述 假设 ， 可 列 出 控制 方程 ， 初 始 条 件 和 
O 流动 并 与 周围 土壤 进行 换 热 ， 根据 测量 得 到 埋 管 进 NF 
o) s zs (ra Sr < oo, T » 0) 
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N 其 中 ,7 为 m 处 的 温 2 穷 举 搜索 法 
E, €; To HERI PEE, °C; ra 为 外 和 孔 
半径 ,mi q 为 埋 管 每 延 米 换 热 量 ，W,m- 1 
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全 GEHERBULHGR, JmK; k 为 土壤 导 邹 系 数 ，。 内 的 不 同 涩 值 组 全 带 入 到 式 (3) 中 ， 得 到 循环 水 温 
全 WE)! 7 为 时 间 , s. 度 随 时 间 变化 的 计算 值 ， 然 后 利用 式 (4) 计算 循环 
= 循环 水 平均 温度 Ti 与 钻 孔 壁 温 Tv. 的 关系 式 ， 温度 计算 值 与 测量 值 之 间 的 方差 和 ， 从 中 寻找 方 
C -— 。\ 差 和 最 小 值 fs， 此 时 对 应 的 土壤 热 物性 参数 值 即 
f^ fiw 一 qiio QR ^p E. EEEE yh ER 

X 认为 是 所 需 测量 的 结果 . 
则 循环 水 进出 口 平 均 温 度 为 I， SN 穷 举 搜索 法 将 热 物性 参数 的 可 能 取 值 全 部 考虑 
og; ner) 在 内 ， 并 逐一 进行 目标 函数 方差 和 的 计算 ， 比 较 得 
Toca t gigei (e OG) — 到 方差 和 最 小 值 。 穷 举 搜索 法 与 优化 方法 相 比 ， 它 


其 中 , Ro 为 钻 筷 内 热 阻 , m KW ,Ei(X) 为 指数 积 
分 函数 。 

参数 估计 法 目标 函数 为 地 下 埋 管 与 周围 土壤 传 
热 模型 计算 得 到 的 循环 水 平均 温度 与 实际 测量 值 之 
间 的 方差 和 M: 


f= DC cali = Toxp, i) (4) 


其 中 , Tou; W i 时刻 由 传 热 模型 计算 得 到 的 地 埋 管 
内 循环 水 平均 温度 ,“C; Texp 为 第 d 时 刻 热 响应 实验 
测 得 的 循环 水 平均 温度 ,“C; N 为 数据 测量 组 数 ， 


可 得 到 目标 函数 最 小 值 ， 而 不 是 某 局 部 区 域 的 极 小 
值 ， 因 而 可 以 得 到 最 接近 实际 的 热 物性 参数 值 ， 穷 
举 搜索 法 不 存在 对 算法 初始 搜索 值 变 化 的 影响 与 依 
赖 。 由 于 穷 举 搜索 法 在 全 局 范围 内 进行 计算 分 析 与 
比较 , 所 以 运算 时 间 较 优化 方法 长 ， 


3 验证 与 分 析 
为 了 验证 穷 举 搜索 法 计算 结果 的 准确 性 ,预先 


设 定 钻 孔 半径 ra, FEK MAE q、 土 壤 导 热 系数 
k, AFLAR Ro、 土 壤 容 积 比 热 p.c. 等 参数 值 ， 
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并 用 以 计算 循环 水 温度 值 作为 模拟 的 测量 结果 ， 然 
后 利用 穷 举 搜索 法 确定 土壤 热 物 性 参数 ， 将 设 定 值 
与 运算 值 进行 对 比分 析 。 采用 两 组 参数 设 定 值 ( 见 
R 1) 分 别 进行 验证 .考虑 到 可 能 的 土壤 热 物性 参 
数值 范围 ， 本 文选 择 的 各 参数 搜索 范围 UU 为 : 0.5 
Wm li.K i<k<5.0 Wm -K^!; 0.01 mK:.W-! < 
Ro < 0.3 mK-W-!, 1.0x10 J.m- 3.K-1l< pscs < 
6.0x105 Jm ?-K-!, 利用 穷 举 搜索 法 的 计算 结果 见 
K2, 

作为 对 比 ， 同 时 采用 优化 方法 寻找 目标 函数 
最 优 解 ， 寻 优 方法 采用 目前 广泛 使 用 的 遗传 算 
法 223), 遗传 算法 通过 “种 群 ” 的 变异 、 交 又 、 选 


择 ， 进 行 非 确定 性 方式 寻 优 ， 最 终 选 定 全 局 范围 内 
所 涉及 “种 群 ” 的 某 一 最 优点 . 为 保证 遗传 算法 计算 
结果 的 准确 性 , 同时 考虑 节省 运行 时 间 , 设置 遗传 算 
法 “遗传 代数 ”为 5000 代 , 每 代 “ 种 群 ”包含 50 个 
个 体 , 优化 方法 计算 得 到 的 结果 见 表 3. 为 进一步 验 
证 两 种 计算 方法 结果 的 可 信和 度 ， 将 模拟 测试 温度 曲 
线 与 两 种 方法 得 到 的 计算 温度 曲线 进行 了 对 比 ， 见 
图 1~2。 

由 表 2 和 表 3 可 知 ， 穷 举 搜索 法 得 到 每 组 方差 
和 为 0, 即 可 以 寻找 到 目标 函数 最 小 值 ,其 对 应 的 各 
参数 值 与 设 定 值 之 间 的 相对 误差 也 均 为 0， 说 明 利 
用 穷 举 搜索 法 可 以 准确 地 找到 热 物性 参数 值 。 而 在 


表 1 设 定 参 数 


Table 1 Preset param 4o 


k/W-m-1-K-1 pscs /106J:m -3-K 71 
第 一 组 0.075 55 1.80 Ke 0.100 2.30 
第 二 组 0.075 75 2.60 0.080 3.50 
表 2 穷 到 的 土壤 热 物性 参数 


B 计算 参数 值 六 下 3、 1. 
相对 误差 /%% 2 0 
i923 2.60 


TEXT 0 0 0 
& 
表 3 优化 方法 得 到 的 土壤 热 物 性 ; 


Table 3 Soil thermal parameters obtained by«optimization method 


设 定 参 数组 别 计算 项 目 
B 计算 参数 值 1.80 V 0413 3.11 
Som FIXES 25/96 0 e 13.00 35.22 PPS 
第 二 组 计算 参数 值 2.61 xS 0.080 3.48 DUE 


相对 误差 /多 0.3 0 0.57 


一 一 模拟 测试 温度 曲线 
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(a) 穷 举 搜索 法 
(a) Exhaustive search method 


0 20 40 60 80 
Es RT/h 

(b) 优化 方法 

(b) Optimization method 


图 1 循环 水 温度 随时 间 变 化 曲线 (第 一 组 ) 


Fig. 1 Circulating water temperature variation vs. time (the first group) 
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一 一 模拟 测试 温度 曲线 


* 穷 举 搜索 法 计算 温度 


0 20 40 60 80 
时 间 /h 
(a) 穷 举 搜索 法 


(a) Exhaustive search method 
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温度 /C 


0 20 40 60 80 
时 间 /h 
(b) 优化 方法 
(b) Optimization method 


图 2 循环 水 温度 随时 间 变 化 曲线 (第 二 组 ) 


Fig. 2 Circulating water temperature variation vs. time (the second group) 


模拟 验证 的 两 组 数据 中 ,优化 方法 得 到 的 结果 与 设 

定 值 之 间 存 在 一 定 误差 ， 其 中 得 到 的 导热 系数 相对 
TO 误差 较 小 , 钻 孔 内 热 阻 相对 误差 较 大 , 可 达到 13%, 
二” ”而 容积 比 热 的 相对 误差 最 大 , 最 高 可 超过 35%, 显然 
优化 方法 得 到 的 导热 系数 准确 度 较 高 ， 但 得 到 的 钻 
孔 内 热 阻 和 容积 比热容 则 可 信 度 相对 较 差 ， 且 优 从 
e eer d 


方差 和 ， 这 表明 优化 方法 寻找 到 了 局 部 J 热 


O 物性 参数 极 小 值 而 不 是 全 局 内 的 最 小 余 |) 关 图 1~2 
N — Bog. PERO USE PRUEBA COM UL 
feu 


C i s 线 和 测试 温 
OQ  。 度 曲线 拟 合 良好 并 不 一 定 表明 计算 结果 一 定 可 靠 . 
> 穷 举 搜索 法 可 以 得 到 确定 的 热 物性 参数 值 ， 且 
C 0 方差 和 最 小 , 计算 温度 曲线 对 模拟 测试 温度 曲线 拟 
C 。 合 非常 好 ， 显 然 穷 举 搜 索 法 得 到 的 热 物性 参数 值 与 
性 “优化 方法 得 到 的 结果 相 比 具有 更 高 的 可 信和 度 . 优化 

方法 寻找 到 的 通常 是 局 部 区 域 极 小 值 而 不 一 定 是 全 


所 以 优化 方法 并 未 实现 真正 意义 上 的 全 局 范围 搜索 
寻 优 . 


答应 用 


AS 
X 利用 穷 举 搜索 法 分 别 对 山东 省 烟台 市 (工程 1) 


和 河北 省 清河 县 (工程 2) 的 两 个 工程 项 目 进行 了 地 
下 土壤 热 物 性 参数 测量 计算 . 现场 测试 得 到 的 进出 
口 循环 水 温度 随时 间 的 变化 见 图 3~4 中 的 测试 温度 
曲线 。 另外， 本 文 还 采用 优化 方法 对 地 下 土壤 热 物 
性 参数 进行 了 计算 , 两 种 方法 得 到 的 结果 见 表 4, 两 
种 方法 利用 土壤 热 物 性 参数 计算 得 到 的 循环 水 温度 
见 图 3~4 中 的 计算 温度 曲线 . 

穷 举 搜索 法 加 关 化 方法 均 得 到 了 较为 接近 的 导 
热 系 数值 , 供 了 内 热 阻 信 存 在 一 定 偏差 ， 而 容积 比 
REZEK. 虽然 两 种 方法 得 到 的 循环 水 计算 温 
度 曲 线 守 实际 测试 温度 曲线 均 拟 合 较 好 ， 但 两 个 工 

穷 举 搜索 法 得 到 方差 和 最 小 值 比 优化 方法 得 到 


局 区 域 最 小 值 ， 其 原因 在 于 遗传 算法 变异 、 交 多 .的 方差 和 极 小 值 均 要 小 ， 显 然 穷 举 搜索 法 得 到 的 热 


过 程 中 没有 得 到 的 “个 体 ”并 未 参与 计算 、 虑 交 


时 间 人 /h 

(a) 穷 举 搜索 法 
(a) Exhaustive search method 
图 3 循环 水 温度 随时 间 变 化 曲线 (工程 1) 


Fig. 3 Curve of circulating water temperature changes with time (The first project) 


物性 参数 值 具有 更 高 的 可 信人 性， 


温度 /C 


时 间 /h 
(b) 优化 方法 
(b) Optimization method 
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穷 举 搜索 法 计算 温度 
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(b) Optimization method 


图 4 循环 水 温度 随时 间 变 化 曲线 (工程 2) 


Fig. 4 Circulating water temperature vs. time (The second project) 


表 4 穷 举 搜索 法 与 优化 方法 得 到 的 土壤 热 物 性 参数 


Table 4 Soil thermal properties obtained by exhaustive se 


计算 方法 


穷 举 搜 索 法 1.73 
super 优化 方法 1.79 
穷 举 搜索 法 1.35 
J 优化 方法 1.34 
5 结 论 


i 


本 文 提出 运用 穷 举 搜索 法 寻找 


进而 确定 土壤 热 物 性 参数 值 ,@ 方法 计算 结 
果 进 行 对 比分 析 , 通过 上 述 耸 


1) 通过 验证 分 析 ， 穷 举 搜索 法 可 以 得 到 目标 函 
数 最 小 值 ， 进 而 确定 准确 的 热 物性 参数 值 ， 其 得 到 
的 计算 温度 曲线 对 模拟 测试 温度 曲线 拟 合 较 好 ， 计 
算 结果 具有 较 高 的 可 信 性 . 优化 方法 得 到 的 土壤 导 
热 系数 较为 接近 预 设 值 ， 但 得 到 的 铬 孔 内 热 阻 及 容 
积 比 热 与 预 设 值 存在 一 定 误差 ， 
穷 举 搜索 法 得 到 的 方差 和 大 . 

人 
的 地 下 土壤 热 物性 参数 进行 计算 , 结果 表明 穷 举 搜 
索 法 较 优 化 方法 得 到 的 地 下 土壤 热 物 性 参数 值 更 为 
ufa. 
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